Negative heat capacity and non-extensive kinetic theory 
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Abstract 

The negative nature of the heat capacity Cy of thermodynamically isolated self-gravitating 
systems is rediscussed in the framework of a non-extensive kinetic theory. It is found that the 
dependence of Cy on the non-extensive parameter q gives rise to a negative branch with the 
critical value corresponding to g = 5/3 (Cy — oo). 
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Consider a simple analogy between ideal gas and self-gravitating system. The tempera- 
ture of this system is defined by 
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-mv'^ = -keT, (1) 

where m is the mass of a particle (e.g., a star) and ks is the Boltzmann's constant. If such a 
system is composed by particles, its total kinetic energy is = ^NksT which, according 
to the virial theorem, is equal to minus its total energy, i.e., E = —K. Therefore, the heat 
capacity of the system is 

which is clearly a negative quantity. In other words, this means that by losing energy 
the system automatically grows hotter, or still, that a self-gravitating system can never 
be in a true thermodynamical equilibrium state (see and references therein). In actual 
fact, it is well known that when discussed in the context of the microcanonical ensemble a 
self-gravitanting system interacting via long range forces exhibits a negative heat capacity 
P]. The negativeness of the heat capacity has been largely studied in many astrophysical 
phenomena. For example, Lynden-Bell & Wood used such a concept to explain the 
gravothermal catastrophe for globular clusters while Beckenstein and Hawking used 
it in the context of black hole thermodynamics. 

On the other hand, a considerable theoretical effort has been concentrated on the study 
of the statistical description of a large variety of physical systems, resulting in the extension 
of the Boltzmann-Gibbs' statistical mechanics. Particular examples are systems endowed 
with long duration memory, anomalous diffusion, turbulence in pure-electron plasma, self- 
gravitating systems or, more generally, systems endowed with long range interactions. In 

n 

order to deal with such problems, Tsallis !6| proposed the following generalization of the 
Boltzmann-Gibbs (BG) entropy formula 

where Pi is the probability of the i*^ microstate and g is a parameter quantifying the 
degree of non-extensivity. In the limit q 1 the celebrated BG extensive formula, 
Si = —kB^^iPi^^Piy is readily recovered as a particular case. 

In the astrophysical context, the first applications of this non-extensive g-statistics studied 



stellar polytropes [3] and the peculiar velocity function of galaxy clusters ^ 



The Jeans 



gravitational instability criterion for a colissionless system was also recently rediscussed 
in the context of this enlarged formalism j9|. In this latter study, it was shown that such 
systems present instability even for wavelenghts of the disturbance smaller than the standard 
critical Jeans value. Recently, Taruya & Sakagami 11 1 discussed a stability problem 



in the microcanonical ensemble by extending Padmanabhan's [3] classical analysis of the 
Antonov instability to the case of polytropic distributions functions. More recently, Chavanis 
[l^ lisl . Il^ considered the Tsallis' entropy as a particular iJ-function corresponding to 
isothermal stellar systems and stellar polytropes. In this approach, the maximization of the 
i?-function at fixed mass and energy reveals naturally a thermodynamical analogy with the 
study of the dynamical stability in collionless stellar systems 24 1 . 

In this Letter, we discuss the negative nature of the specific capacity Cy of self-gravitating 
systems in the framework of a kinetic theory resulting from this non-extensive formalism. 
To do so, we consider self-gravitating systems that can be treated as isolated systems in a 
first approximation, e.g., globular cluster [l^ or elliptical galaxies Q|. We derive a new 
analytical expression for this quantity which gives rise to an entire negative branch. We 
show that for this kind of system the value g = 5/3 is an upper limit for the non-extensive 
parameter. In connection with stellar polytropes, the value of the non-extensive parameter 
g = 5/3 corresponding to a polytropic index n = 3orn=— 1, depending on the expression 
between n and q. 

It has been shown that the kinetic foundations of the Tsallis' themostatistic lead to a 
velocity distribution for free particles given by 
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Mv) = B, 



. mv 



(4) 



2kBT 

In this velocity g-distribution, the g-parameter quantifies the nonadditivity property of the 
associated gas entropy whose the main effect at the level of the distribuition is to replace the 
standard Gaussian form by a power law. Mathematically, this result follows directly from 
the well known identity limw->o(l +wx)^/" = exp(x) jisj. The quantity Bg is a g-dependent 
normalization constant whose expression for g < 1 and g > 1 are respectively given by 

where T is the temperature and n is the particle number density. This non-extensive velocity 



distribution can be derived at least from two different methods: (i) tfirougfi a simple non- 
extensive generalization of the Maxwell ansatz, f{v) ^ f{vx)f{vy)f{vz), which follows from 
the introduction of statistical correlations between the components of the velocities ^| 
and (ii) by using a more rigorous treatment based on the non-extensive formulation of the 
Boltzmann if-theorem, which requires g > jlS^ . 

In order to investigate the g-dependence of the specific capacity, we first consider a cloud of 
ideal gas within the nonrelativistic gravitational context, which is analog to a self-gravitating 
collissionless gas. The kinetic energy of this system is simply given hj K = |m < v"^ >, 
with the statistical content included in the average square velocity of the particles < v"^ >. 
. n this way, the non-extensivity is introduced through a new derivation of expectation value 
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.2 ^ J ~V,n •> 



where = ( m{i-g) ) ^ thermal cutoff on the maximum value allowed for the velocity 
of the particles (g < 1), whereas for the power law g-distribution without cutoff (g > 1) 
Vm ^ oo. This g-expectation value can be easily evaluated for g 7^ 1, and the result 
corresponds to a g-parameterized family of the square velocity of the particles < f ^ >q given 
by 

<v'>,= -^^, for g< 5/3, (7) 
5 — 3g m 

with < >= being the classical value easily obtained in the limit g ^ 1 (It is worth 
mentioning that the virial theorem is not modified in this non-extensive formalism l21]). 

Now, by substituting the non-extensive g-expectation value given by Eq. (0) into the 
heat capacity relation Eq. (j2I) one finds 

Cv = --^Nks, (8) 
5 — 3g 

which clearly reduces to the Maxwellian limit Cy = —\NkB for g = 1. To better understand 
the physical implications of the above expression, in Fig. 1 we plot the quantity Cy/Nks 
as a function of the non-extensive parameter g. The graph clearly shows that the specific 
capacity is a negative quantity for the interval g < 5/3. For g = 5/3 (critical value) the 
specific capacity diverges {Cy —00). From this analysis, it is possible to conclude that 
this particular g-dependence of Cy leads to the classical gravothermal catastrophe {Cy < 0) 
only for power laws with non-extensive parameter in the range < g < 5/3. 




FIG. 1: The quantity Cv/Nks for a self-gravitating collisionless gas obeying the non-extensive 
Tsahis' g-statistic is shown as a function of the non-extensive parameter q. We see that there 
is a negative branch for the interval < g < 5/3. For (7 = 5/3 the specific capacity diverges 
{Cv ^ -oo). 

At this point we emphasize that the true nature of the non-extensive g-parameter appear- 
ing in the Tsalhs' statistical framework remains as a completely open question at present. 
In this way, it is important to show empirical manifestations of these non-extensive effects. 
In this Letter we have discussed a clear signature of non-extensive effects on some kind of 
self-gravitating systems. In particular, we have shown that the g-dependence of the specific 
capacity for these systems is physically consistent only for power-laws lying in the interval 
< q < 5/3. This particular range of values for the non-extensive parameter q is less 
restrictive than that one obtained by Abe js^. We suspect that such a difference occurs 
because the Abe's approach is addressed in the context of the canonical ensemble and, as 
has been shown elsewhere, in this ensemble self-gravitating systems do not present negative 
specific heat 12>]. We also note that, when identified with the polytropic index = | + ^ 
7|, the critical value q = 5/3 is identical to the one that gives rises to the gravitational 
stability/instability transition, i.e., n = 3_ Q| (or n = — 1 if n = ^ + ^|)- In a couple 



of papers by Taruya & Sakagami 

iiii,Q it was shown that a thermodynamic instability 



appears for values of n > 5 in systems confined in adiabatic walls (microcanonical ensem- 
ble) and for values of n > 3 in systems surrounded by a thermal bath. In fact, a more 
consistent study of the gravothermal catastrophe must consider a self-gravitating ideal gas 
contained in a spherical container. Wc shall return to this discussion in a more complete 
analysis involving general non-extensive kinetic theories which will appear in a forthcoming 
communication. 
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